QUANTUM PHASE TRANSITIONS 



D.BELITZ 

Department of Physics and Materials Science Institute 
University of Oregon 
Eugene, OR 97403 

AND 

T.R.KIRKPATRICK 

Institute for Physical Science and Technology 
and Department of Physics 

University of Maryland, College Park, MD 2074-2 



Abstract. These are notes for lectures delivered at the NATO ASI on 
Dynamics in Leiden, The Netherlands, in July 1998. The main concepts 
relating to quantum phase transitions are explained, using the paramagnet- 
to-ferromagnet transition of itinerant electrons as the primary example. 
Some aspects of metal-insulator transitions are also briefly discussed. The 
exposition is strictly pedagogical in nature, with no ambitions with respect 
to completeness or going into technical details. The goal of the lectures is 
to provide a bridge between textbooks on classical critical phenomena and 
the current literature on quantum phase transitions. Some familiarity with 
the concepts of classical phase transitions is helpful, but not absolutely 
necessary. 



Contribution to: Dynamics: Models and Kinetic Methods for Non- 
equilibrium Many-Body Systems 
John Karkheck (editor) 

To be published by Kluwer academic publishers b.v. 



Lectures delivered by D. Belitz 



2 



D.BELITZ AND T.R.KIRKPATRICK 



1. Introduction to Continuous Phase Transitions 

In these lectures we discuss so-called quantum phase transitions or zero- 
temperature phase transitions. Our motivation is the fact that these transi- 
tions provide a field where the long-range correlations in quantum systems 
that were introduced in the preceding lectures [1] have some spectacular 
consequences. Since quantum phase transitions have many features in com- 
mon with ordinary classical or thermal phase transitions, we will first review 
the basic concepts of both, and then define and discuss the particular fea- 
tures of quantum transitions, before turning to some specific examples. We 
would like to warn the reader that it is not possible to do this subject jus- 
tice within the time and space constraints of these lectures. Our aim can 
therefore only be to give some examples that highlight a few of the interest- 
ing features of the field. Thorough treatments of the phenomenology and 
theory of classical phase transitions can be found in Refs. [2], and aspects 
of quantum phase transitions are reviewed in Refs. [3, 4]. 

1.1. REVIEW OF BASIC CONCEPTS 
1.1.1. Continuous phase transitions 

Continuous phase transitions 1 are characterized by the system under con- 
sideration undergoing a transition from a symmetric or disordered state, 
which incorporates some symmetry of the Hamiltonian, to a broken-sym- 
metry or ordered state, which does not have that symmtry, although the 
Hamiltonian still possesses it. A good example is given by a Heisenberg 
ferromagnet: The relevant symmetry of the Hamiltonian is the rotational 
symmetry in spin space, and the disordered and ordered states are repre- 
sented by the paramagnetic and the ferromagnetic states, respectively. In 
the former, there is no preferred direction for the magnetic moments of 
the spins. The net magnetization is therefore zero, and the state has the 
same spin rotational symmetry as the Hamiltonian. That is, if we rotate 
the coordinate system in spin space, then the system will look the same. In 
the latter, on the other hand, there is a preferred direction for the spins, 
and this is the direction which the overall magnetization will point in. The 
state therefore does no longer respect the spin rotational symmetry, and 
we say that the symmetry is spontaneously broken ('spontaneously', in or- 
der to distinguish this phenomenon from the explicit or external breaking 
that occurs if we apply an external magnetic field). In an ideal system the 
preferred direction is completely arbitrary, and in real systems it is deter- 
mined by very small external fields, like e.g. the earth's magnetic field, or 

x We will deal only with continuous transitions, and will simply refer to them as 'phase 
transitions' or 'critical points'. 
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by a small breaking of the symmetry in the system's Hamiltonian as pro- 
vided by, e.g., the ionic lattice strucuture. As a result, in an ideal system 
the magnetization in the ordered state is not only non-zero, but it is also 
non-unique. Such a thermodynamic quantitity that is zero in the disordered 
phase, and non-zero and non-unique in the ordered phase, is called an order 
parameter, a concept first introduced by Landau. The external field B that 
couples to the order parameter, in our example the magnetic field, is called 
the field conjugate to the order parameter. 

As we approach the phase transition by changing some parameter, sev- 
eral remarkable phenomena are observed. For instance, the correlations of 
the order parameter become long-ranged. Let M be the order parameter, 
and let us consider the spatial correlations of M (in the case of a magnet, 
they can be measured by neutron scattering): 

<M(x)M(y))=/(|x-y|/0 . (1) 

Everywhere in parameter space except at the critical point, the function 
/ decays in space on a length scale £ called the correlation length. 2 £ 
diverges as the critical point is approached, usually like a power law that 
is characterized by the correlation length critical exponent, u: 

, (2) 

where t denotes some dimensionless distance in parameter space from the 
critical point. For instance, if the transition occurs at a non-zero critical 
temperature T c , we can use t = \T — T c \/T c . At criticality, i.e. at t = 0, 
the correlation length diverges, which indicates that the order parameter 
correlations decay only like a power law, which has no intrinsic scale. 3 It 
is customary to denote this power law by an exponent r/: 

(M(x)M(y)) t=0 cx|x-y|- d - 2 ^ , (3) 

where d denotes the spatial dimensionality of the system. 

Apart from these long-range correlations in space, there are similar ef- 
fects in the temporal behavior of the system. Let us denote the equilibration 
time, i.e. the time scale for the system to return to equilibrium after it has 
been disturbed, by r c . This equilibration time diverges as criticality is ap- 
proached, and it does so as a power of the correlation length, with the 
power law characterized by an exponent z: 

T C ^e , (4) 

2 For most phase transitions, the function / in Eq. (1) decays exponentially for large 
arguments: In f(x — ► oo) = —x + 0(lna:). However, in general the functional form can be 
more complicated. 

3 Mathematically, a power law is a homogeneous function, while an exponential, for 
instance, is not. 
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The inverse of r c defines a critical frequency scale lo c that goes to zero as 
criticality is approached, a phenomenon called critical slowing down: 

oj c (t -» 0) oc 1/t c -> . (5) 

The exponents u, r], and z that we have defined so far are examples 
of critical exponents, which characterize the power law behavior of various 
observables upon approach to the critical point. Three other important crit- 
ical exponents are (3, which describes the vanishing of the order parameter, 

M(t -»• 0) oc iP , (6) 

7, which describes the diverges of the order parameter susceptibility, which 
in our example is the magnetic susceptibility \ = M/B, 

X(*-0)oct-T , (7) 

and the critical exponent 5, which describes the dependence of the order 
parameter on its conjugate field at criticality, 

M(t = 0, B -» 0) oc B l/S , (8) 

Notice that a diverging susceptibility implies that there cannot be a region 
where the order parameter depends linearly on the conjugate field. 

The set of critical exponents characterizes the critical behavior, and it 
turns out that all of the critical exponents are not independent. Rather, 
many of them are related to one another by scaling laws or exponent rela- 
tions. Furthermore, it turns out that the complete set of critical exponents 
is the same for whole classes of phase transitions. For instance, the criti- 
cal exponents for the ferromagnetic transitions in iron and nickel are the 
same, even though these three metals have very different bandstructures, 
and their Curie temperatures are very different. Perhaps even more sur- 
prisingly, the critical exponents observed at the critical point of water are 
the same as those at the liquid-gas critical points of the "quantum liquids" 
He3 and He4. Phase transitions that share the same critical behavior, as 
expressed by the critical exponents, are said to belong to the same uni- 
versality class, and the existence of such universality classes is referred to 
as universality or universal behavior. It turns out that universality classes 
are determined by basic symmetries of the underlying Hamiltonian, and by 
the spatial dimensionality of the system. The fact that the critical behavior 
is independent of the microscopic details of the Hamiltonian is due to the 
diverging correlation length: Close to a critical point, the system performs 
an average over all length scales that are smaller than the (very large) cor- 
relation length. This observation also gives a clue as to how to theoretically 
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deal with critical phenomena: In order to correctly describe the universal 
critical behavior it should be sufficient to work with an effective theory 
that keeps explicitly only the asymptotic long-wavelength or large-distance 
behavior of the original Hamiltonian. By contrast, observables that vary 
from system to system even within a given unversality class, like e.g. the 
critical temperature, are called non-universal properties, and they are not 
accessible by means of such effective theories. 

The theoretical understanding of these, and other, properties of phase 
transitions took place over a period of 100 years, starting with van der 
Waals, continuing with Landau, and ending with Wilson. Wilson's renor- 
malization group, the essence of which is a set of scale transformations, 
provided us with a deep understanding of the origin of universality and the 
critical power laws. More recently, it has also turned out to be a very useful 
general tool for studying the statistical mechanics of many-body systems 
in general, not only close to some critical point. Here it is not our goal to 
describe the renormalization group, for this purpose we refer the reader to 
the many excellent expositions in the literature [5]. Rather, we will take 
the critical phenomenology for granted, and on its basis discuss the role of 
quantum mechanics in the context of phase transitions. The only exception 
is Sec. 2.3, where we assume some familiarity with basic renormalization 
group arguments in order to derive some results. As an explicit example, we 
will stick to ferromagnets, except in Sec. 3 which deals with metal-insulator 
transitions. 

1.2. CLASSICAL VERSUS QUANTUM PHASE TRANSITIONS 

A fundamental question is the following: To what extent is quantum me- 
chanics necessary in order to understand the critical phenomena we have 
sketched in the previous subsection, and to what extent will classical physics 
suffice? One can get the correct answer by means of the following very sim- 
ple considerations (which can be elaborated upon if that is considered desir- 
able). Generally speaking, quantum mechanics is important whenever the 
temperature becomes lower than some characteristic energy of the system 
under consideration. For instance, in an atom that characteristic energy is 
the Rydberg energy. In our case, we have seen that there is a character- 
istic frequency, namely uj c . Let us assume that the corresponding energy 
scale, huj c , is the smallest relevant energy scale. Since u c (t — > 0) — > 0, this 
is a reasonable guess close to the transition. It then follows that quantum 
mechanics will be important whenever the temperature T obeys 

k B T < huj c . (9) 

Usually we think of phase transitions as taking place at some non-zero criti- 
cal temperature T c , e.g. at the Curie temperature in our magnetic example. 
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It then follows that sufficiently close to the transition, namely for 



we expect quantum mechanics not to be important for describing the sys- 
tem's behavior. Here To is some microscopic temperature scale, e.g., the 
Fermi temperature in a metallic ferromagnet. It follows that for any phase 
transition that takes place at a non-zero critical temperature T c > 0, the 
critical behavior asymptotically close to the transition can be described 
entirely by classical physics. This conclusion survives more rigorous argu- 
ments. These phase transitions are called classical or thermal transitions. 
What drives the correlation length to infinity are thermal fluctuation, which 
become very large close to criticality. In contrast, we might think of a tran- 
sition that occurs at zero temperature, and that is triggered by varying 
some non-thermal parameter, e.g. the system's composition. For instance, 
if n is the concentration of some ingredient, and n c is the critical concen- 
tration, we can choose t = \n c — n\/n c as our dimensionless distance from 
criticality. Then it obviously follows that quantum mechanics will be im- 
portant for describing the critical behavior. These transitions are called 
quantum or zero-temperature phase transitions, and the relevant fluctuata- 
tions are quantum fluctuations or zero-point motion. An example would 
be a magnetic system that is, at T = 0, continuously diluted with some 
non-magnetic material until it undergoes a transition to the paramagnetic 
state. (We will deal with the obviously relevant question of what happens at 
low, but non-zero, temperatures in the following subsection.) Notice that, 
according to this definition, some phase transitions in systems that are 
usually considered quintessentially quantum mechanical, like the supercon- 
ducting transition in mercury at T = 4.2 K, or the A-transition in helium 
at T = 2.17 K, are classical transitions. Indeed, in both cases the critical 
behavior (but not the physics that triggers the transition, or the properties 
of either phase) can be understood entirely by means of classical physics: 
Both transitions are in the universality class of a classical 3-d xy-model. 

1.2.1. Classical phase transitions 

Let us now first consider some additional properties of classical phase tran- 
sitions. Within classical statistical mechanics, consider a Hamiltonian 




(10) 



H(p,q)=K(p) + U(q) 



(11) 
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where p and q are the generalized momenta and positions, and K and U 
are the kinetic and potential energy, respectively. 4 The partition function 

Z = J dpdq e- H ^ T = J dp e~ K ^ T J dq e~ U ^ T , (12) 

then factorizes into a piece that depends only on K and one that depends 
only on U. As a result, one can study the system's static properties in- 
dependently from its dynamical ones. In particular, the dynamical critical 
exponent z is independent from all of the other critical exponents, and 
the static critical behavior can be studied, following Landau, by means of 
an effective functional of a time- independent order parameter. One often 
expresses this by saying that 'statics and dynamics decouple'. 

Close to the critical point, the free energy density, / = (— 1/V) k^T In Z, 
obeys a generalized homogeneity law, 

f(t,B,...) = b- d f(tb l ^,Bb^,...) . (13) 

Here V is the system volume, and t and B are the dimensionless distance 
from the critical point and the external field conjugate to the order param- 
eter, respectively, as before, b is an arbitrary positive real number called a 
scale parameter, and Eq. (13) holds for all b > 0. v is the correlation length 
critical exponent, and xb > is a critical exponent that is related to 5 
by xb = dS/(l + 5). Since all thermodynamic quantities can be obtained 
from the free energy, Eq. (13) provides us with homogeneity laws for all 
of them. For instance, by differentiating / with respect to B we obtain a 
homogeneity law for the order parameter density, M = df/dB, 

M(t, B) = b XB - d M(t b l/v , Bb XB ) . (14) 

Since b is arbitrary, we can in particular put b = B~ l / XB . At t = we then 
recover the above relation between xb and 5. By the same method, other 
exponent relations can be obtained, and it turns out that of all the static 
critical exponents, only two are independent, e.g. v and xb- Another useful 
substitution is to set b = t~ u oc £. This makes it obvious that letting b — ► oo 
is tantamount to approaching criticality. In this context, a remark about the 
suppressed variables in Eq. (13), which we denoted by is in order. 

They all enter Eq. (13) analogously to t and B, but the exponents that 
characterize their 'scaled' entries on the right-hand side (i.e. the analogs of 
\jv and xb), turn out to be negative. As a result, these entries go to zero as 
one approaches criticality, and are called 'irrelevant variables' or 'irrelevant 

4 We exclude from our considerations systems of charged particles in magnetic fields, 
and other cases of velocity dependent potentials. 



8 



D.BELITZ AND T.R.KIRKPATRICK 



operators'. If the observable under consideration is a regular function of a 
particular variable for small values of the argument, then it follows that 
this variable becomes unimportant as we approach criticality and does not 
influence the critical behavior. It thus really becomes 'irrelevant' in the 
ordinary sense of the word. This is not the case, however, if the observable 
is a singular function of some argument for small values of the argument. 
In this case the irrelevant variable influences the critical behavior after all, 
and one speaks of a 'dangerous irrelevant variable'. Dangerous irrelevant 
variables can cause substantial complications in the technical analysis of 
scaling near critical points. 

The homogeneity or scaling law, Eq. (13), was historically first postu- 
lated phenomenologically, and it turned out that all observed properties of 
phase transitions followed from it if appropriate values for v and xb were 
used, depending on the universality class under consideration. It was the 
triumph of Wilson's renormalization [6] group that it allowed a derivation 
of the homogeneity law from first principles. This derivation is highly tech- 
nical, and we cannot go into it. Instead, we refer the reader to the extensive 
literature on this subject [5]. 

1.2.2. Quantum phase transitions 

Let us now turn to the case of quantum phase transitions, i.e. transitions 
that occur at T = and are triggered by some non-thermal control pa- 
rameter. As we have seen above, in this case quantum mechanics is always 
important, and we need to employ quantum statistical mechanics in order 
to calculate the partition function and the free energy. Let the Hamilton 
operator of the system be H(a^,a), with and a a collection of creation 
and annihilation operators. In the usual imaginary time formalism, the time 
evolution of a is given by d(r) = e Hr ae~ Hr , with r the imaginary time 
variable. A general theorem [7] then tells us that the partition function for 
any quantum many-body system can be written as a function integral of 
the form 

Z = J D[i>^} e 5 ^ . (15) 

Here ?/> and ip are space and imaginary time dependent fields that are 
isomorphic to the sets of creation and annihilation operators in a second 
quantization formulation of the problem. Their nature depends on whether 
the quantum particles are fermions or bosons: For the latter, the fields 
are classical or bosonic (i.e., they commute), while for the former, they 
are fermionic (i.e., they anticommute) . D[ip,ip] is an appropriate integra- 
tion measure defined with respect to these fields [8]. The action Sfy,^] is 
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uniquely determined by the Hamilton operator, and is given by 

d 



-I 



r ri/k B T 
S[4>,^} = dx dr ■0( x J r 

l/k B T 



V>(x,r) 



drff(^(x,T)i(x,r)) . (16) 



Here H as a function of z/> and ■0 has the same functional form as H as a 
function of d^ and a. 

Since the Hamiltonian taken at some imaginary time does not commute 
with the Hamiltonian taken at another imaginary time, we see that for 
quantum systems the statics and the dynamics are intrinsically coupled and 
need to be treated together and simultaneously. All phenomena that need 
to be described by means of quantum statistical mechanics, and quantum 
phase transitions in particular, therefore automatically fall under the title 
of this NATO AST It further follows that for quantum phase transitions, 
in contrast to classical ones, there are three independent critical exponents, 
and the dynamical exponent z needs to be determined together with the 
static ones. 

The homogeneity law for the free energy density, Eq. (13), can now 
easily be generalized to the quantum case. From Eqs. (2) and (4,5) we see 
that, as a function of i, frequencies scale like t uz . Furthermore, in the imag- 
inary time formalism, temperature and Matsubara frequencies are directly 
proportional to one another, and it is therefore plausible that temperature 
and frequency will scale in the same way. We thus add T as an argument 
to our free energy density and acknowledge the explicit T in the definition 
of / to obtain 

f(t,T,B,...) = b-^ f{tb 1/l/ ,Tb z ,Bb XB ,...) . (17) 

Comparing Eqs. (13) and (17) we see that a quantum phase transition in d 
spatial dimensions resembles the corresponding classical transition in d e g = 
d + z spatial dimensions! This is also plausible from the point of view of the 
Landau functional (see below) , where a spatial integral J <ix in the classical 
case gets replaced by a space-time integral J dxdr in the quantum case. 
Early work on the subject suggested that this observation provides a fast 
and easy solution to the problem of quantum critical behavior. However, 
as we will see, the argument is too superficial to be reliable, and the extent 
to which it holds requires a careful and detailed discussion. 
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Figure 1. Schematic phase diagram showing a paramagnetic (PM) and a ferromagnetic 
(FM) phase. The path indicated by the dashed line path represents a classical phase 
transition, and the one indicated by the solid line a quantum phase transition. 



1.3. AN EXAMPLE: THE PARAMAGNET-TO-FERROMAGNET 
TRANSITION 

Let us illustrate the concepts introduced above by means of a concrete 
example. For definiteness, we consider a metallic or itinerant ferromagnet. 5 

1.3.1. The phase diagram 

Figure 1 shows a schematic phase diagram in the T-J plane, with J the 
strength of the exchance coupling that is responsible for ferromagnetism. 
The coexistence curve separates the paramagnetic phase at large T and 
small J from the ferromagnetic one at small T and large J. For a given J, 
there is a critical temperature, the Curie temperature T c , where the phase 
transition occurs. This is the usual situation: A particular material has a 
given value of J, and the classical transition is triggered by lowering the 
temperature through T c . Alternatively, however, we can image changing J 
at zero temperature (e.g. by alloying the magnet with some non-magnetic 
material). Then we will encounter the paramagnet-to- ferromagnet transi- 
tion at the critical value J c . This is the quantum phase transition we are 
interested in. Since we have seen that the quantum transition is, loosely 
speaking, related to the classical one in a different spatial dimensionality, 

5 For the purposes of this subsection we might as well consider localized spins, but the 
theory discussed in Sec. 2 below applies to itinerant magnets only. 
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Figure 2. Schematic phase diagram as in Fig. 1 showing the vicinity of the quantum 
critical point (J = J c , T = 0). Indicated are the critical region, and the regions dominated 
by the classical and quantum mechanical critical behavior, respectively. A measurement 
along the path shown will observe the crossover from quantum critical behavior away 
from the transition to classical critical behavior asymptotically close to it. 



and since we know that changing the dimensionality usually means chang- 
ing the universality class, we expect the critical behavior at this quantum 
critical point to be different from the one observed at any other point on 
the coexistence curve. 

This brings us to the question of how continuity is guaranteed when one 
moves along the coexistence curve. The answer, which was found by Suzuki 
[9], also explains why the behavior at the T = critical point is relevant 
for observations at small but non-zero temperatures. Consider Fig. 2, which 
shows an enlarged section of the phase diagram near the quantum criti- 
cal point. The critical region, i.e. the region in parameter space where the 
critical power laws can be observed, is bounded by the two dashed lines. 
It then turns out that the critical region is divided into two subregions, 
denoted by 'QM' and 'classical' in Fig. 2, in which the observed critical 
behavior is predominantly quantum mechanical and classical, respectively. 
The division between these two regions (shown as a dotted line in Fig. 2) 
is not sharp (and neither is the boundary of the critical region) , but rather 
a smooth crossover from predominantly quantum mechanical to predom- 
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inantly classical critical behavior is observed as J is increased at low but 
non-zero T. The figure illustrates that the asymptotic critical behavior, very 
close to the transition, is classical for all non-zero T, but that for small T 
there nevertheless is a sizeable region where quantum critical behavior is 
observable. It also makes it clear that the abovementioned continuity is 
realized by means of a distribution limit. 

1.3.2. Classic results for itinerant ferromagnets 

The subject of quantum magnetic phase transitions was pioneered by Hertz 
[10], who built on earlier work by Suzuki [9] and others [11]. Hertz's main 
result can be recovered by combining the above observation that the quan- 
tum phase transition should correspond to the classical one in d e s = d + z 
dimensions, with the fact that for classical magnets the upper critical di- 
mension d£, i.e. the dimensionality above which the critical behavior is 
mean-field like, is d£ = 4. It then follows that for the quantum transition, 
d+ = 4 — z. Hertz further found that z = 3 for clean itinerant ferromagnets, 
and z = 4 for disordered ones (we will reproduce his results in Sec. 2 below). 
It then seemed to follow that d£ = 1 for clean magnets, and d£ = for 
disordered ones, so that the quantum critical behavior would be mean-field 
like in all physical dimensions (d = 2, 3). 

These classic results, which seemed to solve the problem (albeit the an- 
swer was not very interesting from a theoretical point of view), cannot be 
correct, however, since they violate other, very general results that were ob- 
tained around the same time. Curiously, this contradiction went unnoticed 
for more than 20 years. 

In 1974, Harris [12] argued that at any critical point in a system that 
contains quenched disorder, the correlation length exponent must obey the 
inequality 

v > 2/d , (18) 

in order for the critical behavior to be stable with respect to the disorder. 
6 Harris' physical arguments were later augmented by Chayes et al. [13], 
who proved a rigorous mathematical theorem to the same effect. Now the 
mean-field value of v is vmf = 1/2, which violates Harris' criterion for all 
d < 4. Hence Hertz's results for disordered systems cannot be correct. It 
was also pointed out by Sachdev [14] that Hertz's results for clean systems 
in d < 1 (an academic, but nonetheless interesting case) were at odds with 
some general scaling arguments. 

It finally turned out that the classic results for both the clean and 
the disordered case are indeed incorrect, for rather subtle and interesting 

6 This is not how Harris formulated his criterion, but it can be brought into this form 
by using suitable exponent relations. 
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reasons. It also was shown that the theory can be salvaged with relatively 
little effort, provided that one acknowledges the long-range correlations in 
itinerant electron systems that were explained in the preceding lectures [1]. 
This theory was developed in a series of papers [15], the main arguments 
and results of which we reproduce in the next section. 

2. Quantum Critical Behavior of Itinerant Ferromagnets 

We now sketch the derivation of an effective field theory for itinerant fer- 
romagnets, starting from a microscopic Hamiltonian. As we will see, the 
derivation makes massive use of our knowledge of the properties of non- 
magnetic electron systems. It is this knowledge that makes the task of de- 
termining the correct quantum critical behavior much easier than it would 
have been in the 1970s. 

2.1. LANDAU-GINZBURG- WILSON THEORY 

Let us consider a microscopic model for an itinerant ferromagnet. What 
we have in mind is to develop an effective theory from such a miscroscopic 
model that will be capable of correctly describing the quantum critical 
behavior without being more detailed than necessary in other respects. To 
this end we follow Landau's strategy [17] of formulating a theory in terms 
of an appropriate order parameter, and integrating out all other degrees of 
freedom. For a ferromagnet, the relevant observable is the spin density 

n s (x,r) = 4>(x,T)aip(x,T) , (19) 

where the vector a denotes the Pauli matrices. The relevant electron- 
electron interaction is the spin-triplet interaction between the spin density 
fluctuations, 

r ri/ksT 

S'int = J j dx j dT n s( x > r ) ' n s(x,r) . (20) 
The complete action is 

S = So + S'int • (21) 

Here So comprises all parts of the action other than Si n t as defined in Eq. 
(20): It contains a term describing free electrons (or lattice electrons 7 ) as 
well as electron-electron interactions in all channels other than the spin- 
triplet one. In particular, So contains the direct or Coulomb interaction 
between the electronic charge densities. We will refer to the fictitious system 

7 For our purposes, which are aimed at long-wavelength phenomena, it is irrelevant 
whether one considers a lattice model or a continuum one. 
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that is described by the action So as the reference ensemble, and it will play 
an important role in what follows. 

Our goal is to calculate the partition function 



D S {i(>,ip] g J / dx n s (x)-n s (x) ^2) 



Here and it what follows we adopt a four- vector notation x = (x, r) for 
space-time integrals. We now concentrate on the term Si nt , and decouple it 
by means of a Gaussian or Hubbard-Stratonovich transformation [18]. The 
latter consists of introducing a classical (i.e. commuting) auxiliary field 
M{x) that couples linearly to the spin density, and rewriting e Sint as a 
Gaussian integral, 

Z = J D[ip,ip] e So ^ J d[M] e~ J f dx M2(x)+2J f dx Mix)ns(x) . (23) 

From now on, for simplicity we ignore the vector nature of the spin den- 
sity. Taking it into account just complicates the notation without lead- 
ing to qualitatively important effects. Next, we interchange the M and 
■0-integrations, and formally carry out the latter, 

Z = J D[M] e~ J I dx M2{x) Z [M] , (24) 

where we recognize 

Z [M] = J D$,l/>] e So$M+2Jf** M(x)n s (x) ^ ( 25 ) 

as the partition function of the reference ensemble in an external 'magnetic 
field' that is proportional to M(x). From the linear coupling between M{x) 
and the spin density it is clear that the expectation values of these two 
quantities are proportional to one another. We thus identify M{x) as the 
order parameter field whose expectation value is the magnetization. It is 
customary to define a Landau-Ginzburg- Wilson or LGW functional $[M], 

$[M] = J J dx M 2 (x) - In Z [M] , (26) 

where 

Z [M] = {e 2J I dx M (*W*))„ . (27) 

Here (. . .)o denotes an average with respect to the reference ensemble ac- 
tion So, and the partition function can be expressed in terms of the LGW 
functional as 

Z = J D[M] e~*W . (28) 
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The above derivation of an order parameter theory for quantum ferro- 
magets is due to Hertz [10]. The salient point is that the LGW functional 
is given in terms of the free energy of the reference ensemble in a (space 
and time dependent) 'external field' that is given by the order parameter 
field. A more technical way to say this is that the LGW is given by the 
generating functional for connected spin density correlation functions in 
the reference ensemble, viz. InZo. Notice that all of the above have been 
exact manipulations. The strategy underlying this exact rewriting of the 
partition function is to express the free energy of the full system, which 
undergoes a phase transition, in terms of that of the reference ensemble, 
which does not. 8 The price one pays is that one needs to know the free en- 
ergy of the latter in the presence of an arbitrary space and time dependent 
magnetic field. As we will see, however, enough is known about correlated 
electrons in magnetic fields to successfully implement this strategy. 

2.2. THE LANDAU EXPANSION 

Follwing Landau, the next step is to expand the LGW functional in pow- 
ers of M. Remembering that the reference ensemble has no spontaneous 
magnetization, and hence (n s )o = 0, we have 

Z = l + 2J 2 J dxdy M(x)M(y)(n s (x)n s (y)) + ... , (29) 

where the terms not shown explicitly are of 0(M 4 ) or higher. Taking a 
Fourier transform, we thus obtain for the LGW functional 

q ui 

+ ]T u 4 ({q,^})M 4 (q,u,) + ... , (30) 

{q,"d 

where q = |q|, and 

Xs(q,u) = (n s (q,w)n s (-q, , (31) 

is the connected two-point spin-density correlation function of the reference 
ensemble, i.e., its spin susceptibility, 

U4 = (n s n s n s n s )o , (32) 

8 Remember that the reference ensemble is missing the spin-triplet part of the electron- 
electron interaction that is responsible for triggering ferromagnetic phase transitions. 
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is the corresponding connected four-point function, etc. We have used an 
obvious schematic notation for the quartic terms, since we will not study 
them in detail here. 

The salient point of this formal development is that the LGW functional 
is given in terms of the connected spin density correlation functions of the 
reference ensemble. The reference ensemble, however, is a Fermi liquid, or 
its generalization to the case of quenched disorder, and so its correlation 
functions are known! Indeed, the preceding lectures [1] has explained in 
some detail what they are, and we can now draw on that knowledge. For 
the sake of definiteness we will discuss only disordered systems, 9 where the 
effects we want to demonstrate are most pronounced. Qualitatively similar, 
albeit weaker, phenomena are present in clean systems as well, as has been 
discusssed in the original literature [15]. The appropriate limit to study 
the correlation functions is that of small wavenumbers and frequencies, 
q,u> — ► 0, with u <C q 10 in suitable units. 11 As we have seen in the 
preceding lectures [1, 20] the spin susceptibility of a disordered Fermi liquid 
in this limit reads 

Xs(q,u J )= X ° s (q) [l-H/Dq 2 + ...] , (33) 
and the long-wavelength expansion of the static spin susceptibility Xs reads 
X° s (q) = const. -q d ' 2 -q 2 + {q 2 ) , (34) 

where o(q 2 ) denotes terms that are smaller than q 2 , and we have omit- 
ted positive prefactors of all terms in the expansion, as they will not be 
important for what follows. The LGW functional now takes the form 



$ [ M ] = 77E T E M (q> w ) \t + q d ~ 2 +q 2 + MAZ 2 1 M(-q, 



V Q 



-L0 



+0(M 4 ) . (35) 



The four-point correlation function U4 also contains nonanalyticities, which 
show up the in the quartic term in <£ [M] , and the same is true for all higher 

9 Strictly speaking, we should have used a replicated theory to deal with the quenched 
disorder. In the interest of keeping our pedagogical discussion simple we have suppressed 
this technical point. For an introductory discussion of the replica trick, see Ref. [19]. 

10 Otherwise one does not reach criticality, which can be seen as follows. Since the 
magnetization is conserved, ordering on a length scale L requires some spin density to 
be transported over that length, which takes a time t oc L 2 / D, with D the spin diffusion 
coefficient. Now look at the system at a momentum scale q or a length scale L oc 1/q < £, 
with £ the coherence length. Because of the time it takes the system to order on that scale, 
the condition for criticality is L 2 < Min(Di, £ 2 ). In particular, one must have L 2 < Dt, 
or lu oc 1/t < Dq 2 . 

11 For instance, one can measure q in units of the Fermi wavenumber fcp, and ui in units 
of Dkp, with D the spin diffusion coefficient of the reference ensemble. 
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order terms. While a careful study of these terms is necessary for a complete 
treatment of the problem, we suppress them here for brevity and simplicity, 
and refer the interested reader to the original literature where they have 
been discussed in detail [15]. 

Before we analyze the LGW functional, Eq. (35), in the next subsec- 
tion, two remarks are in order. First, if we had used a reference ensemble of 
non- interacting electrons, rather than our more realistic one that includes 
electron-electron interactions, then we would have missed the non-analytic 
terms proportional to q d ~ 2 in Eqs. (34) and (35). As a consequence, we 
would have recovered Hertz's results which, as we have seen earlier, can- 
not be correct. While it seems at this point as if a realistic choice of the 
reference ensemble were crucial, this is disturbing from some fundamental 
theoretical points of view. Indeed, a careful investigation reveals that the 
precise choice of the reference ensemble is not important. We will come back 
to this point in Sec. 2.4 below. Second, the physical origins of the q d ~ 2 term 
are long-range spatial correlations in the reference ensemble of interacting 
electrons that underlies our effective LGW action, as has been explained in 
the preceding lectures [1]. These long-range spatial correlations are in turn 
a consequence of the coupling of statics and dynamics in quantum systems. 
We now see what has happened: In addition to the order parameter fluctu- 
ations that develop a long range near criticality, our electron systems also 
contains long-range correlated degrees of freedom (viz. the 'diffusons' of the 
preceding lectures) that have nothing to do with phase transition physics, 
and that are present even far away from the transition. These degrees of 
freedom have been integrated out in our derivation of the LGW functional, 
which by definition is a functional of the order parameter field only. The 
inevitable consequence of this integrating out of slow modes is nonanalytic- 
ities in the resulting LGW functional. Such nonanlyticities violate the very 
spirit of the LGW concept, and to perform explicit calculations for the 
non-local field theory given by Eq. (35) would be very difficult indeed. An 
obvious way to avoid these problems would be to not integrate out the dif- 
fusons, but rather derive a generalized LGW theory in terms of all the soft 
modes in the systems. However, as we will see in the next subsection, for 
the purpose of determining the critical behavior the non-local theory can 
be handled, and the present route is the fastest one to answer the questions 
we have asked. 

2.3. RENORMALIZATION GROUP ANALYSIS 

To finish our treatment of itinerant ferromagnets, we now analyze the LGW 
functional, Eq. (35), by means of power counting or a tree- level renormal- 
ization group (RG) analysis. Space and time constraints do not allow us to 
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explain this technique here. Readers not familiar with it can find excellent 
and very accessible treatments in Refs. [2]. 

Let us restrict ourselves to spatial dimensions d < 4, where the q d ~ 2 
term dominates over the analytic q 2 term, and let us look for a Gaussian 
RG fixed point where neither the q d ~ 2 term nor the \u>\/q 2 term in Eq. 
(35) are renormalized. At such a fixed point, the critical order parameter 
correlation function will behave like 

(MM) t=0 <xl/q d - 2 ^q 2 ^ , (36) 

where the last relation reflects the definition of the critical exponent rj, see 
Eq. (3). We thus have 

77 = 4 - d . (37) 

Furthermore, at such a Gaussian fixed point the frequency clearly scales 
with the wavenumber like u> ~ q d = q z , 12 which yield a dynamical critical 
exponent 

z = d . (38) 

Finally, the exponent 7 for the Gaussian fixed point can also just be read 
off the Gaussian action: At zero frequency and wavenumber, we have 

(MM)^=o«l/t = r , (39) 

so that we have 

7 = 1 . (40) 

In order to determine the correlation length exponent v, we define the scale 
dimension [q] of a wavenumber q to be [q] = 1. The requirement that the 
q d ~ 2 term in the LGW functional be dimensionless then leads to a scale 
dimension of the order parameter field of [M] = — (d — 2)/2. Since the t 
term must also be dimensionless, this yields 

t ~ q d ~ 2 ~ r ( d - 2 ) = C 1/v , (41) 

from which we read off v as 

v = l/(d - 2) . (42) 

All of these results obviously hold only for 2 < d < 4. For d < 2 the 
electrons in the reference ensemble become localized (see the preceding 
lectures) and our theoretical framework breaks down, while for d > 4 the 

12 We use oc for 'proportional to', and ~ for 'scales like' or 'has the same scale dimension 
as'. 
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q 2 term dominates over the q d ~ 2 term, and we recover Hertz's mean-field 
exponents, i.e. r\ = 0, z = 4, 7 = 1, and v = 1/2. 

A quick check shows that our result fulfils the Harris criterion, v > 2/d 
(see Sec. 1.3.2 above) for all values of d. While this is encouraging, it is 
of course only a necessary criterion for our results representing the correct 
quantum critical behavior, not a sufficient one. To establish the latter, 
one needs to consider the higher order terms in the LGW functional, and 
establish that our Gaussian fixed point is stable. At tree level, this has been 
done in the original literature [15], and the result was that the fixed point is 
indeed stable. The analysis of the quartic term in $ in particular also yields 
the equation of state, and thus the critical behavior of the order parameter 
itself, i.e. the critical exponents (3 and 8. The result is 

P = 2/(d-2) , 6 = d/2 , (43) 

for 2 < d < 6, while for all d > 6 these exponents have their mean-field 
values (3 = 1/2 and <5 = 3, respectively. The analysis of the higher order 
terms thus establishes d = 6 as another upper critical dimensionality, in 
addition to d = 4. In d = 4 and in d = 6 logarithmic corrections to scaling 
occur, as is usually the case at an upper critical dimensionality. 

Finally, we mention again that an analogous analysis of clean systems 
yields qualitatively very similar results. Essentially, the exponent d — 2 in 
the non-analytic terms gets replaced by d— 1, which leads to a single upper 
critical dimension d~£ = 3. Contrary to the disordered case, the critical 
behavior in the most interesting dimension d = 3 is therefore mean-field 
like with logarithmic corrections to scaling [15]. 

2.4. FINAL REMARKS 

We finish this section with a few additional remarks. 13 As we have sketched 
in the preceding subsection, and as has been more carefully established in 
the original literature, we have managed to determine the critical behavior 
at the quantum ferromagnetic transition exactly, yet this behavior is not 
mean-field like. This is surprising, as usually non-mean field like critical 
behavior cannot be obtained exactly, save for a very few models. To find 
out what has enabled us to do so, let us look again at our LGW functional, 
Eq. (35). One way to state the effect of the non-critical soft modes that 
have led to the q d ~ 2 term is to say that they have established an effective 
long-range interaction between the order parameter fluctuations. Indeed, 
a Fourier transform of the non-analytic term yields an interaction that 
falls off like l/r 2<yd ~ l \ It is well known that such long-range interactions 

13 Like the preceding subsection, some of these remarks require familiarity with renor- 
malization group techniques. 
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stabilize Gaussian critical behavior that is not mean-field like [21]. What is 
remarkable here is that this long-range interaction is not put in by hand, 
but rather is generated by the system itself via the non-critical slow modes. 

These considerations raise the question of whether the phenomenon 
discussed here is germane to quantum phase transitions. In principle, it is 
not. Whenever there are slow modes in addition to the order parameter 
fluctuations that couple to the latter, one will obtain nonanalyticities in 
the LGW functional, and resulting unusual critical behavior, irrespective of 
whether one deals with a classical or a quantum phase transition. However, 
quantum transitions are much more susceptible to this mechanism, for the 
simple reason that there are many modes that are soft only at T = and 
acquire a mass at non-zero temperature. Our diffusons are a good example: 
As was discussed in the preceding lectures, they are indeed massive at 
T > 0. 

Finally, we come back to a point raised at the end of Sec. 2.2. There 
we mentioned that whether or not one obtains the correct critical behavior 
seems to depend crucially on the choice of the reference ensemble. While 
one might argue that a non-interacting reference ensemble is simple not a 
realistic model, this leads to the following paradox: Suppose we consider a 
model whose action consists only of a free electron part, and the spin-triplet 
interaction S- m t of Eq. (20). Then a decoupling of S- m t, as performed above, 
seems to lead to mean- field critical behavior. However, if we considered some 
fraction of S[ n t a part of So, an d decoupled the rest, then we would have 
an interacting reference ensemble and would obtain the above non-mean 
field like critical behavior. Clearly, both procedures are equally valid and 
should lead to the same result. The resolution of this paradox is as follows. 
In the case of an interacting reference ensemble, already the bare LGW 
functional contains the crucial nonanalyticities. Therefore, a RG analysis 
at tree level is sufficient to obtain the correct critical behavior. In the case 
of a non-interacting reference ensemble, on the other hand, the bare action 
does not contain the crucial non-analytic terms, but they are generated if 
the RG analysis is carried to higher order in the loop expansion. Indeed, an 
inspection of Hertz's model shows that the q d ~ 2 term is indeed generated by 
the RG, starting at one-loop order. Since the generated term is relevant with 
respect to the mean-field fixed point, it invalidates the zero-loop analysis. 
These observations serve as a reminder of a fact that is well-known in 
principle, but occasionally forgotten: Any RG analysis to a given order in 
a loop expansion gives the correct answer only if no relevant new terms in 
the action are generated at higher order. Unfortunately, proving that no 
such terms are generated at any order in the loop expansion amounts to 
proving that the theory is renormalizable, a task that is very difficult and 
has been done only for a few select models. 
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3. The Anderson-Mott Transition 

At the end of these lectures, we would like to briefly touch upon some very 
different quantum phase transitions, namely the types of metal-insulator 
transition known as Anderson, and Anderson-Mott transitions, respectively. 
Our motivation for doing so is chiefly to dispel any possible misconception 
that the concepts developed so far apply only to quantum phase transitions 
that are magnetic in nature, or that are represented as the T = end point 
of a line of classical phase transitions. Indeed, the Anderson-Mott transition 
is neither. We would also like to illustrate the point that the scaling ideas 
that have been so successful in connection with phase transitions, classical 
and quantum mechanical, may be applied to transport coefficients as well as 
to thermodynamic quantities [22]. Apart from making these two point, the 
time and space we devote to this subject are grossly inadequate. Extensive 
reviews can be found in Refs. [23, 3, 24, 25]. 

3.1. THE SIMPLEST METAL-INSULATOR TRANSITION: THE CLASSICAL 
LORENTZ MODEL 

The simplest metal-insulator transition occurs in a classical model, viz. the 
classical Lorentz gas. It is well known that with increasing scatterer density 
n, the diffusivity D of the moving particle decreases (see the preceding 
lectures), until if finally reaches zero at a critical scatterer density n c . This 
behavior can be seen in the numerical data shown in Fig. 3. With t = 
\n — n c \/n c the dimensionless distance from the critical point, one has 

D(t ^0)(xt s , (44) 

with s a critical exponent. The reason for the localization of the diffusing 
particle is that with increasing scatterer density it is increasingly likely to 
get trapped in cages of scatterers from which it cannot escape. 

3.2. ELECTRONS: THE ANDERSON TRANSITION, AND THE MOTT 
TRANSITION 

Let us now add quantum mechanics to these considerations, and consider 
the quantum Lorentz model that was discussed in detail in the preceding 
lectures. Quantum mechanics has two competing effects on the transport 
properties: On the one hand, it should increase the diffusivity, since it allows 
the particle to tunnel out of cages it would be trapped in classically. On the 
other hand, it leads to the quantum interference or weak-localization effects 
that enhance the backscattering amplitude. As we saw in the preceding 
lectures, the latter effect wins. The quantum Lorentz model therefore has a 
stronger tendency to localize the particle than the classical one, to the point 
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Figure 3. Numerical simulation data for the diffusion coefficient D versus the scatterer 
density n of a 2-d classical Lorentz model. Full symbols represent data by Bruin [30], 
and open ones data by Alder and Alley [31], and by Alley [32]. D is normalized by its 
Boltzmann value I)' ', and n by its critical value n c . From Ref. [23]. 



that in d = 2 the particle is localized even for arbitrarily small scatterer 
density. In d > 2, however, there is a metal-insulator transition at a finite 
value of n, which is called an Anderson transition. 

The Anderson transition is a strange phase transition from a statistical 
mechanics point of view, as it has no simple order parameter, and no up- 
per critical dimensionality. However, it has a lower critical dimensionality, 

= 2, as for d < 2 there is no transition. This has been exploited for 
studies of the Anderson transition in e = d — 2 expansions. More generally, 
the dynamical conductivity a (or, alternatively, the diffusion coefficient D) 
obeys a generalized homogeneity law [26] 

tr(t,u) =b-( d ~ 2) a(tb 1/u ,ujb z ) . (45) 

Putting u = 0, and b = t~ u , we obtain from Eq. (45) the behavior of the 
static conductivity, 

a(t,0)(xt s , (46) 
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with a conductivity exponent 

s = u(d-2)>2(d-2)/d . (47) 

The equality in Eq. (47) is known as Wegner's scaling law, and the inequal- 
ity results from the Harris criterion, Eq. (18). Due to the poor convergence 
properties of the 2 + e expansions, no reliable theoretical values for s or 
v are available. Numerical calculations in d = 3 yield values for v in the 
range 1.3 - 1.5 [25]. 

The quantum Lorentz model does not contain any electron-electron in- 
teraction, and the Anderson transition is therefore entirely driven by disor- 
der. The opposite case, namely a metal-insulator transition that is entirely 
driven by interactions, with no quenched disorder present, was proposed by 
Mott to explain why certain materials with one electron per unit cell, e.g. 
NiO, are insulators. Mott's original idea hinged on the long-range nature 
of the Coulomb interaction, and in this case the metal-insulator transition 
comes about by means of a breakdown of screening and is of first order [27]. 
A similar, albeit continuous, transition is believed to occur in a model with 
a short-ranged electron-electron interaction known as the Hubbard model. 
This Mott-Hubbard transition is still not well understood in d = 3, al- 
though much progress has been made recently on high-dimensional models 
[28]. 

3.3. THE ANDERSON-MOTT TRANSITION 

An Anderson-Mott transition results if one combines the driving forces for 
the Anderson and Mott transitions and considers the case of interacting 
electrons in the presence of disorder. This is a problem of great experimen- 
tal interest. For instance, the metal-insulator transition that is observed 
in doped semiconductors as a function of the dopant concentration is be- 
lieved to be of that type. Fig. 4 shows the infrared absorption spectrum 
of phosphorus-doped silicon. For very low dopant concentrations one sees 
a hydrogen-like spectrum that is produced by isolated phosphorus atoms. 
With increasing dopant concentration these 'atoms' start to overlap, which 
leads to a broadening of the spectral features. At the highest concentra- 
tion shown in the figure, the spectrum is smooth, but it still represents an 
insulator (no absorption at zero energy) . With further increasing donor con- 
centration, the system undergoes a quantum phase transition to a metal, 
as can be seen in Fig. 5, where both the static conductivity and the dielec- 
tric susceptibility are plotted versus the phosphorus concentration. This 
system, and some other doped semiconductors, are the experimentally best 
studied examples of metal- insulator transitions. Notice that this is a quan- 
tum phase transition that has no classical counterpart, as there can be no 
true insulator at any non-zero temperature. 
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Figure 4- Infrared absorption coefficient a for three different donor concentrations nu 
in Si:P as measured by Thomas at al. [33]. The critical concentration in this system is 
n c « 3.7 x fO 18 cm" 3 . From Ref. [23]. 



Much theoretical effort has been devoted to the Anderson-Mott transi- 
tion. These approaches fall into two distinct classes. The first one contains 
e-expansions about the lower critical dimension d~ = 2 [23]. These theo- 
ries lead again to Wegner scaling, Eq. (45), as in the case of the Anderson 
transition, and Wegner 's scaling law, Eq. (47), still holds. This leads to 
s = v > 2/3 in d = 3, which contradicts the experimental result that the 
value of s is close to 1/2. 14 

The second class is formed by an order parameter theory [3] that uses 
the fact that the density of states at the Fermi level may be considered 
an order parameter for the Anderson-Mott transition (in contrast to the 
case of the Anderson transition, where the density of states is uncritical.) 
This theory has established d£ = 6 as the upper critical dimension for the 

14 All experimentalist do not agree with this result. Ref. [29] has reported s > 1 for 
Si:P. 
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Figure 5. Divergence of the dielectric susceptibility \ (full symbols), and vanishing of 
the static conductivity a (open symbols), both extrapolated to zero temperature, at the 
metal-insulator transition in Si:P as measured by Rosenbaum et al. [34]. From Ref. [23]. 



Anderson-Mott transition (again in contrast to the Anderson transition, 
where = oo), and the critical behavior for d > 6 is exactly known 
and mean-field like. Certain technical analogies between this theory and 
theories for magnets in random magnetic fields have led to the suggestion 
that in d < 6, and in particular in d = 3, the Anderson-Mott transition has 
aspects that are reminiscient of a glass transition, with exponential rather 
than power-law critical behavior for many observables. A scaling theory has 
been developed for this unorthodox critical behavior [3], which is at least 
consistent with existing experimental results, but so far no microscopic 
theory exists. 
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